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We compute chiral logarithms in the presence of "taste" symmetry breaking of staggered fermions. The 
lagrangian of Lee and Sharpe is generalized and then used to calculate the logs in tt and K masses. We correct an 
error in Ref. [|^]; the issue turns out to have implications for the comparison with simulations, even at tree level. 
MILC data with three light dynamical flavors can be well fit by our formulas. However, two new chiral parameters, 
which describe 0{a^) hairpin diagrams for taste-nonsinglet mesons, enter in the fits. To obtain precise results for 
the physical 0{p*) coefficients, these new parameters will need to be bounded, at least roughly. 



It has become clear that simulation at rather 
light quark mass is crucial for accurate determi- 
nation of physical parameters, e.g., heavy-light 
decay constants [||. Since staggered fermions 
provide the fastest known method of simulating 
low-mass quarks, the systematic effects associ- 
ated with this fermion choice should be studied. 

We adopt the following nomenclature: A single 
staggered field describes 4 equivalent "tastes" of 
quarks in the continuum limit; taste symmetry is 
broken at O(a^). We use the word "flavor" for 
different staggered fields; lattice flavor symmetry 
is exact for equal masses. MILC's improved stag- 
gered ("a^-tad") simulations |^ use three fields 
("u, d, s") and reduce the tastes to one per fla- 
vor by taking v^Det. We call these "2 -I- 1" flavor 
simulations since we take m„ = md = me. 

MILC data for vs. quark mass show clear 
deviations from linearity, as expected from chi- 
ral logarithms. However, the detailed behavior at 
a « 0.13 fm does not agree with continuum chi- 
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ral perturbation theory (xPT) Q|. Figure ^ shows 
a fit to the continuum forms for / (mi -I- 7712), 
where X = tt or and toi, m2 are quark masses. 
(We define "pions" to have toi = m2 and "kaons" 
to have m2 ~ m^^^^^ and mi 7^ m2.) The fit is 
very poor, with a confidence level of 5 x 10"^. 

To get good fits, one needs to include the 0{a^) 
taste-breaking effects into the chiral calculations, 
i.e., we need "staggered chiral perturbation the- 
ory" (SxPT). Such calculations start with the 
chiral lagrangian of Lee and Sharpe |^ for 1 stag- 
gered flavor (4 tastes). One defines 

E(a;) = exp(i0//) ; = ^5^5 + <^mCp + ' ' ' , (1) 

where </> is a 4 x 4 matrix of pseudoscalar mesons 
of various tastes, and ^5, . . . are taste matrices. 
The Lee- Sharpe lagrangian is then 

£ - :^tr(9^Sa^I]t)-^^tr(I] + S]t) + a2v (2) 
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The taste-breaking potential V has various op- 
erators with explicit taste matrices; —V = 
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Figure 1. ^^/(mi +m2) vs. mi +m2, in units 
of the potential scale ri, for 2 + 1 flavor lattices at 
a ~ 0.13 fm. The flt is to the continuum chiral log 
forms. The upper branch for (mi +m2)ri ~ 0.2 uses 
the kaon form; the rest of the curve is the pion form. 

SiLi CiOi. For example, 

O2 = i[tr(ES])-tr(e5Se5S)+h.c.] (3) 
O5 = ^[tr(^,Se.S^)-tr(^,^5Se5e.St)] . (4) 

To get NLO chiral results relevant to 2 + 1 sim- 
ulations, we follow a three step procedure: 

1. Generalize the Lee-Sharpe lagrangian to the 3- 
flavor case, i.e., include three lattice staggered fields 
{u,d,s), each with four tastes, and masses rriu ~ 
■nid = me 7^ nis. This is an "8 + 4" theory. 

2. Compute the desired quantities (here, mj and 
m|f ) at one loop in the 8 + 4 case. 

3. Adjust the 8 + 4 answer "by hand" to correspond 
to the 2 + 1 case of interest by identifying the chi- 
ral contributions that correspond to n virtual quark 
loops and dividing them by 4". 

Step 1 turns out to be non-trivial: Rcf. Q cm- 
ployed Fierz transformations to simplify the fla- 
vor structure of the generalized potential V. How- 
ever, Ref. already used the Fierz freedom in 
deriving V. The net result is that the generalized 
operators O2 and O5 are incorrect in Ref. 

There are actually two operators that corre- 
spond to each of O2 and O5: 

021 = i[Tr(e^S)Tr(C^E)+h.c.] 

022 = i[Tr(e^^5S)Tr(^5CMS) + h-C.] 



051 = ^Tr(e^E)Tr(e^I]t) 

052 = iTr(e^C5S)Tr(e5epSt) . (5) 

In the absence of flavor indices, the combinations 
O21 + O22 and — O51 + O52 can be Fierzed into 
Lee-Sharpe O2 and O5, respectively; the other 
linear combinations turn into other Oi. 

The new operators have a surprising effect: 
they generate hairpin (disconnected) diagrams for 
taste- nonsinglet (but flavor neutral) mesons. For 
taste vector and taste axial vector mesons, this 
occurs at chiral tree level. The magnitudes of the 
hairpins are 

S'y^^{C21-C5iy, S'^^^{C22-C,2) (6) 

for vector and axial taste, respectively. 

The disconnected hairpin diagrams have not 
been included in any simulations to date of taste- 
nonsinglet mesons. Thus the Lee-Sharpe la- 
grangian does not apply to such simulations: The 
chiral lagrangian requires more than one flavor to 
describe particles, e.g., tt"*" — ud, that have no 
disconnected contributions. 

At 1-loop, we need to iterate three kinds of 
hairpins on flavor-neutral internal lines: the stan- 
dard anomaly hairpin (toq) for taste singlets, and 
the new hairpins for taste vector and axial vector. 
Unlike TOq, Sy cannot be taken to infinity. We 
therefore have to rediagonalize the mass matrix in 
the taste V, A channels. We call, for example, the 
mass eigenstates in the flavor-neutral, taste vec- 
tor channel iTy , rjv and rj'y . Techniques in Ref. ^ 
are very useful for reexpressing the iterated prop- 
agator as a sum of simple poles. 

Adjustment to go from the 8 + 4 to the 2 + 1 
theory is easy. Every hairpin interaction on a 
meson line after the first introduces an additional 
virtual quark loop. Therefore, if -F(A) is the sum 
of hairpin diagrams, just take F(A) 4F(|-) for 
A = S'y, S'^ or 5 (where S = ml/{2ATT^p)). 

Additional diagrams that were identified as va- 
lence diagrams in Ref. turn out to be discon- 
nected and contribute only to flavor-neutral cor- 
relators. 

We deflne (3j^+ to be the chiral log term for the 
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Goldstone K'^ mass {m\~l°°^)'^ / {ii{me + iris)) = 
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l + /3^+/(167r2/2) + 
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Our results are then: 



•, and similarly for /3^4 
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Here 5' is the ss meson. The continuum result 
is in the taste-singlet channel (/); the remainder 



vanishes in the limit 5[, 



0. 
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With 8y as free parameters, these results give 
excellent fits to the MILC data. However, the 
parameters are poorly constrained: 5y and 5^ 
wander off in opposite directions to large values 
(an order of magnitude larger than known taste- 
violating terms in the 7r+ sector). On the other 
hand, if we fix 5y or to be of reasonable magni- 
tude, the fits find reasonable values for the other 
parameter, too, and the confidence levels are still 
excellent. Figure ^ shows the fit with S'j^rf fixed 
to -0.1; 5'yrl is then found to be 0.25(11). Based 
on the measured taste violations, the natural size 
for these parameters is 5'r\^Q.2. 

We conclude that SxPT is necessary to fit cur- 
rent staggered lattice data. The taste-violating 
hairpins introduce two new free parameters {^'y a) 
into the chiral theory; other taste violations, 
which split flavor- nonsinglet mesons, are not free 
parameters in the fits because they are measured 
directly in simulations. Unfortunately, the hair- 
pin parameters are not well constrained, at least 
at present. It does however appear to us that even 
a rough constraint on 5y {e.g., demanding that 
they be no more than 3 times the known taste- 
violating terms) will be sufficient to determine 
analytic coefficients (Gasser-Leutwyler "Li") 
with good precision. 

We may be able to constrain the hairpin param- 
eters in a variety of ways. Additional SxPT cal- 
culations are underway that generalize these 
results to the quenched and partially quenched 




Figure 2. Same as Fig. ^, but using the SxPT forms, 
eq. (0). The parameter S'j^ is held fixed to a reason- 
ably small value in the fit. 

cases, and extend them to pseudoscalar decay 
constants (of both light-light and heavy-light 
mesons). With these calculations in hand, the 
direct constraints from all the fits may prove suf- 
ficient. Another approach is to compute the four- 
uark, taste-violating operators perturbatively 
; the hairpin parameters may then be estimated 
by vacuum saturation or lattice calculation of the 
matrix elements. Finally, a direct lattice evalua- 
tion of the disconnected hairpin graphs may be 
possible. 
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